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ABSTRACT 


Computer  programs  have  been  developed  for  calculation 
of  the  molecular  absorption  and  aerosol  extinction  for 
horizontal  path  propagation  of  radiation  in  the  wavelength 
range  above  0.55  ym.  The  AFCRL  line  parameters  compilation 
is  used  as  input  data  with  local  weather  conditions  for  the 
Monterey  Bay  and  model  of  typical  aerosol  size  distribution 
from  Shettle  and  Fenn.  Wind  speed  and  relative  humidity 
are  included  in  the  form  developed  by  Wells,  Gal  and  Munn. 
Infinite  resolution  computations  have  been  carried  out  for 
1.06,  3.8  and  10.6  ym  and  results  tabulated  for  monthly 
and  yearly  average  weather  conditions  on  Monterey  Bay. 
Bandwidth-averaged  spectral  transmittances  have  been  calcu- 
lated at  0.01  cm  1 and  0.5  cm  ^ resolution  for  the  same 
conditions. 

Comparisons  with  previously  published  data  show  agree- 
ment to  1.8%  with  maritime  model  aerosol  Mie-scattering 
calculations  by  Selby  using  LOWTRAN  III  B [Ref.  12] . 

Analysis  of  sensitivity  to  changes  of  pressure  and 
temperature  for  the  3.8  ym  range  showed  only  small  effects 
of  these  parameters.  However  a change  in  water  vapor 
partial  pressure  results  in  a 0.75%/mbar  change  in  trans- 
mittance in  the  3.8  ym  region. 
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I . INTRODUCTION 


A.  BACKGROUND 


The  light  bringing  us  information  from  a distance, 
illuminating  our  surroundings  and  carrying  energy,  per- 
forms these  functions  only  by  travelling  in  the  atmos- 
phere. Laser  light  ir.  the  form  of  a highly  collimated 
beam  with  very  narrow  frequency  bandwidth  and  high  radiant 
power  is  equally  dependent  on  the  propagation  properties 
of  the  atmosphere.  The  characteristics  of  the  light: 
notably  its  color  (wavelength) , intensity,  polarization 
and  spatial  distribution  are  all  modified  by  interaction 
with  the  constituents  of  the  atmosphere  - mostly  the  absorb- 
ing trace  gases  and  suspended  particles  in  the  size  range 
up  to  the  order  of  hundreds  of  micrometers. 

In  modern  military  technology  a wide  range  of  systems 
are  in  use  or  under  development  for  surveillance,  imaging, 
target  detection,  designation,  identification  and  tracking, 
for  communications,  for  precision  guidance  of  munitions 
and  possibly  even  as  a directed  energy  weapon.  These  oper- 
ate over  the  entire  "optical"  range  of  wavelengths,  from 
the  visible  or  near  ultra-violet  through  the  "far  infra  red"  — 
from  about  0.3  to  20  micrometers.  Operation  of  these  sys- 
tems is  limited  by  the  properties  of  the  atmosphere  as  a 
transmission  medium,  particularly  the  degradation  of 
resolution  by  scattering  from  particulates  and  turbulence. 
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and  loss  of  radiant  intensity  through  absorption  and  scatter- 
ing. Due  to  strong  absorption  by  atmospheric  molecules, 
systems  are  constrained  to  operate  in  limited  wavelength 
ranges  or  "windows",  relatively  free  of  absorption  lines. 

Of  these  windows,  the  radiating  temperatures  of  many 
sources  of  military  interest  and  the  desirability  of  day/ 
night  operating  capability  have  concentrated  attention  on 
the  3-5  urn  and  8-14  um  infra-red  windows. 

In  particular  the  navy  has  a special  interest  in  the 
DF  laser  as  a candidate  moderate  or  high  power  laser  at  a 
wavelength  in  the  3-5  um  range,  unaffected  by  the  high 
humidity  levels  found  in  the  naval  environment.  Of  particu- 
lar interest  is  the  3.8007  P2(8)  transition  line  which 
shows  up  to  96%  transmission  over  a kilometer  path  at  sea 
level. 

The  need  to  predict  the  operational  performance  of 
these  systems  leads  to  a requirement  to  predict  the  spec- 
tral transmittance  of  the  atmosphere  in  the  infra-red,  as 
a function  of  the  meteorological  parameters  which  determine 
the  absorption  and  scattering  levels  — the  pressure,  tem- 
perature, air  composition,  vertical  gradients,  number  den- 
sity composition  and  size  distribution  of  suspended  particles. 

Large  scale  measurement  programs  have  been  undertaken 
in  many  countries  with  the  objective  of  building  a large 
library  of  meteorological  and  transmission  data.  Consider- 
able effort  (notably  at  AFGL)  has  been  devoted  to  the 
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development  of  computer  programs  for  transmittance  calcu- 
lation based  on  weather  data,  and  the  establishment  of  a 
number  of  "standard"  atmospheric  models. 

The  current  modelling  effort  was  undertaken  in  support 
of  the  experimental  program  of  transmittance  measurement 
on  ranges  over  Monterey  Bay  being  carried  out  by  the 
Environmental  Physics  Group  at  the  Naval  Postgraduate 
School.  This  work  involves  measurement  in  all  the  suitable 
atmospheric  windows  from  the  visible  to  14  urn,  with  both 
laser  and  broadband  sources . The  calculational  program 
has  therefore  been  directed  to  computation  of  molecular 
absorption  and  aerosol  extinction  coefficients  for  sea 
level  propagation  under  a range  of  weather  conditions  sta- 
tistically typical  of  Monterey.  This  has  been  done  for  a 
variety  of  wavelengths  appropriate  to  the  measurements, 
and  at  high  resolution  for  use  with  laser  sources,  and 
averaged  over  filter  bandwidths  for  broadband  sources.  The 
Air  Force  Cambridge  Research  Laboratory  Atmospheric  Absorp- 
tion Line  Parameters  Compilation  [Ref.  7]  in  magnetic  tape 
form  is  used  as  input  data  for  the  absorption  calculations. 

A Mie  scattering  program  based  on  the  Van  De  Hulst  approxi- 
mation with  the  model  particle  distribution  of  Shettle  and 
Fenn  [Ref.  8]  has  been  used  for  the  aerosol  scattering. 

Wind  speed  and  relative  humidity  effects  have  been  included 
using  the  drop  size  growth  model  of  Wells,  Gal  and  Munn 
[Ref.  18].  A sensitivity  analysis  has  been  carried  out  with 
respect  to  pressure  and  temperature. 


12 


Basically,  this  thesis  describes  a program  to 

a)  Develop  a computer  code  to  calculate  atmospheric 
transmittance  for  specific  wavelengths  in  the  near  and 
mid  IR. 

b)  Develop  a code  to  give  transmittance  and  molecular 
absorptance  with  bandwidth  appropriate  to  grey-body  sources 
in  the  near  and  mid  IR. 

c)  Apply  these  to  prediction  of  transmittance  over 
marine  optical  paths,  over  Monterey  Bay,  based  on  the 
available  meteorological  data  base  for  the  area,  for  con- 
ditions appropriate  to  the  optical  propagation  experimental 
measurement  program. 

Curves  have  been  prepared  giving  predicted  molecular 
and  aerosol  extinctions  for  statistical  monthly  and  yearly 
average  meteorological  conditions  on  Monterey  Bay. 

Further  sections  of  this  thesis  describe  the  theoreti- 
cal background  and  the  details  of  the  algorithms,  and  provide 
tabular  and  graphical  results  of  the  calculations.  Computer 
program  listings  are  given  in  the  appendices . 

B.  ATMOSPHERIC  PHENOMENA 


Optical  transmission  of  energy  is  usually  described 
by  the  transmittance. 


where  y is  called  the  extinction  coefficient,  and  L repre- 
sents the  length  of  the  path  which  the  light  has  traveled. 
This  expression  follows  from  Beer's  Law,  which  for  linear 
propagation  of  monochromatic  radiation  at  frequency  v has 
the  form 


gj  I(z,v)  = -y(v)  I(z,  y)  (1-2) 

Here  I(z,y)  represents  the  intensity  of  monochromatic 
radiation  at  a distance  z. 

Beer's  Law  in  differential  form  (1-2)  may  also  be 
written  as 


I (L,  y)  = KO,  y)  e_w(v)L  (1-3) 

for  a given  path  length  L.  The  transmittance  x over  a 
path  of  length  L is  then  given  from  the  equation  (1-3) 


x (v) 


I (L, v)  - y ( v ) L 

1(0, v)  e 


which  is  the  same  as  equation  (1-1)  for  monochromatic 
radiation. 

The  extinction  coefficient  y is  generally  the  result 
of  two  processes:  absorption  and  scattering.  Both  mole- 
cules and  particles  suspended  in  the  atmosphere  cause 
absorption  and  scattering.  It  may  be  useful  to  define  the 
extinction  coefficient  y as  a sum  of  four  terms 
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V = k + 0 + k + 0 

m m a a 


(1-4) 


where  k , B , k , 8=  are  molecular  absorption,  molecular 
m m a a r 

scattering,  aerosol  absorption,  aerosol  scattering 
coefficients . 

An  aerosol  extinction  coefficient  y is  sometimes  used 

cl 

for  the  summation  of  aerosol  absorption  and  scattering. 

Similarly  a molecular  extinction  coefficient  y is  used 

m 

for  the  summation  of  molecular  absorption  and  scattering. 

From  equations  (1-1)  and  (1-4)  the  total  transmittance 
through  a path  in  the  atmosphere  may  be  written  as 


x (total)  = t (molecular  absorption)  * 
x (molecular  scattering)  * 
x (aerosol  absorption)  * 
t (aerosol  scattering) 


where  x (molecular  absorption)  may  be  written  in  terms  of 


the  atmospheric  components  as 


t (molecular  absorption)  = t (uniformly  mixed  gases  line  abs . ) * 

t (ozone  line  abs.)  * 
t (water  vapor  line  abs.)  * 
t (water  vapor  continuum  abs . ) * 
x (nitrogen  continuum  abs.) 
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Although  the  molecular  scattering,  aerosol  absorption 
and  aerosol  scattering  coefficients  are  smooth  functions 


of  frequency  of  monochromatic  radiation,  the  molecular 
absorption  coefficient  is  strongly  dependent  on  frequency. 

Molecules  which  cause  absorption  and  scattering  of  a 
light  beam  in  the  atmosphere  may  be  divided  into  two  types: 
uniformly  and  non-uniformly  mixed  gases.  Carbon  dioxide, 
nitrous  dioxide,  carbon  monoxide,  methane,  oxygen  and 
nitrogen  are  assumed  as  uniformly  mixed  gases  for  which 
their  mixing  ratios  by  volume  with  each  other  are  constant, 
independent  of  pressure  and  temperature.  These  mixing  ratios 
are  tabulated  in  table  4.  Water  vapor  and  ozone  are  con- 
sidered as  non-uniformly  mixed  gases,  whose  relative  con- 
centrations must  be  specified. 

Six  different  models  of  the  atmosphere  are  frequently 
used  in  calculations;  these  are  the  tropical,  midlatitude 
summer,  midlatitude  winter,  subarctic  summer,  subarctic 
winter  and  U.S.  Standard  atmosphere  1962  model  described  by 
McClatchey  and  D'agati  [Ref.  4].  Pressure,  temperature, 
density,  water  vapor  and  ozone  concentrations  in  the  air 
for  each  model  at  sea  level  are  tabulated  in  table  5 . 

The  number  densities  of  aerosol  particles  which  also 
cause  absorption  and  scattering  of  a light  beam  are  more 
complicated  to  determine.  It  is  necessary  to  know  the  size 
distribution  of  aerosol  particles  in  addition  to  the  number 
density,  in  order  to  compute  the  aerosol  scattering  and 


absorption  coefficients.  The  size  distribution  function 
is  denoted  by  n(r)  and  given  as 


n (r) 


dN  (r) 
dr 


(1-5) 


where  N(r)  is  called  the  number  density  which  is  the  number 
of  aerosol  particles  with  radius  less  than  r per  unit  volume. 

Several  different  analytical  models  have  been  developed 
for  the  size  distribution.  A "continental"  aerosol  model, 

\ 

which  includes  both  rural  and  urban  aerosol  conditions, 
and  a "maritime"  aerosol  model  are  used  for  the  lower 
atmosphere. 

The  rural  aerosol  model  used  by  Shettle  and  Fenn  is 
assumed  to  be  composed  of  a mixture  of  70%  water  soluble 
substance,  consisting  of  ammonium  and  calcium  sulfate  and 
organic  compounds,  and  30%  dustlike  aerosols  [Ref.  8] . 

In  urban  areas  the  rural  aerosol  background  is  modified 
by  the  addition  of  aerosols  from  combustion  products  and 
industrial  sources.  The  proportions  of  soot-like  aerosols 
and  the  rural  aerosol  type  aerosol  mixture  are  assumed  to 
be  35%  and  65%  respectively  [Ref.  8]. 

The  maritime  aerosol  model  represents  the  aerosol  com- 
position and  size  and  number  distributions  over  the  oceans. 

It  is  significantly  different  from  the  continental  aerosol 
types.  The  maritime  model  is  due  to  salt  particles  caused 
by  the  evaporation  of  the  seaspray  droplets  also  added  to 


the  rural  type  aerosols,  from  which  most  of  the  very  large 
particles  have  been  eliminated. 


The  general  form  of  the  size  distribution  for  these 
models  is  given  by  Shettle  and  Fenn  [Ref.  8]  as 


n (r)  = l (- 


N. 


i=1  (In  10)  rai  (2tt) 


(log  r - log  r . ) 

exp  [ - 5 1 


2ai 


(l-6a) 

where  the  parameters  in  the  equation  above  are  tabulated 
in  table  6. 

The  size  distribution  given  by  equation  (l-6a)  is 
relative  humidity  and  wind  speed  independent.  Since  the 
number  of  oceanic  particles  with  greater  radius  in  the  size 
distribution  increases  with  increasing  relative  humidity  and 
wind  speed,  a modified  maritime  model  has  been  proposed 
for  calculation  of  aerosol  extinction  coefficients  taking 
care  of  the  growth  of  particle  radius  [Ref.  18] . This 
gives 

n (r)  = £r{0.47(£r)~A  + 2 . 3MRa  (p-)  exp  [-8 . 5 (p-)  Y]  } (l-6b) 

where  C is  a normalization  constant,  and  MR  is  the  mixing 
ratio  of  seaspray  aerosols  with  the  rural  aerosols.  Or 
more  clearly,  the  mixing  ratio  of  rural-oceanic  aerosols 
in  the  equation  (l-6b)  is  1:MR. 

Particle  growth  is  included  through  the  growth  factor  F 
which  adds  the  effect  of  relative  humidity  into  the  equation 
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(l-6b) . It  is  given  as 


1 - 0.9  In (1  - 


where  R.H.  is  the  percent  relative  humidity.  The  effect 
of  F is  to  increase  the  drop  radius  to 


rR.H.  " F ro 


where  rQ  is  the  radius  for  zero  percent  relative  humidity. 

F is  normalized  to  80%  relative  humidity,  given  F'.  There- 
fore factor  F'  in  the  equation  (l-6b)  is 


F'  = F(R.H.) 

F ( 80 ) 


Y in  equation  (l-6b)  describes  the  dependence  on  the  wind 
speed.  It  has  the  form 


1.25 


y = 0.384  - 0.00293  V 


where  V is  the  wind  speed  in  meters  per  second,  a in 
equation  (l-6b)  also  represents  the  wind  speed  dependence 
and  it  is  given  as 


a = 250  + 750  v1,16  , when  V < 7 m/sec. 


a = 6900  V 


0.29 


when  V > 7 m/sec, 


Z- 


So  the  inputs  for  equation  (l-6b)  are  wind  speed  V and 
relative  humidity  R.H.  The  normalization  constant  C may 
be  determined  from  the  relation  between  visibility  and 
aerosol  extinction  coefficient  at  0.55  micrometers. 

C.  HF/DF  LASER  AT  THE  NAVAL  POSTGRADUATE  SCHOOL 

Although  this  thesis  will  deal  with  propagation  at 
sea  level  at  any  frequency  of  radiation,  primary  interest 
will  be  the  frequencies  of  the  HF/DF  laser.  These  fre- 
quencies are  tabulated  in  tables  7 and  8 . 

To  be  used  for  the  propagation  measurement  program  at 
the  Naval  Postgraduate  School,  a small  scale  hydrogen 
fluoride  (HF)  or  deuterium  fluoride  (DF)  laser  was  located 
on  the  roof  of  Spanagel  Hall.  The  system  consists  of  a 
vacuum  system,  a power  supply  and  the  laser  head.  The 
laser  is  centered  around  a laser  head  which  contains  dis- 
charge tube,  mixing  chamber  or  cavity  and  resonant  cavity. 
The  system  uses  four  gases:  a)  Sulfur  hexafluoride,  b)  oxy- 
gen, c)  helium,  d)  hydrogen  or  deuterium.  A mixture  of 
sulfur  hexafluoride,  oxygen  and  helium  is  injected  into  the 
discharge  tube  at  the  end  opposite  the  mixing  chamber. 

Then  the  sulfur  hexafluoride  is  dissociated  by  electron 
impact  in  the  discharge  tube  by  making  use  of  an  electric 
discharge.  Oxygen  helps  the  sulfur  hexafluoride  to  dis- 
sociate into  free  fluorine  radicals  by  the  reaction 


SF.  +0-  S0_F_  + 2F 

4 2 2 2 


1 


I 

The  hydrogen  or  deuterium  is  injected  perpendicular  to  the 
flow  in  the  mixing  chamber.  Then  chemical  reaction  takes 
place  in  this  mixing  region  forming  the  vibrationally 
excited  HF  or  DF  molecules. 

The  optical  resonator  of  the  system  consists  of  a 
pair  of  Brewster  windows  and  mirrors.  The  Brewster  windows 
which  linearly  polarize  the  output  are  chosen  from  two  inter- 
changeable pairs,  one  of  BaF2  (0.25-15  vim  transmitting) 
and  the  other  sapphire  (0.14-6.5  ym  transmitting). 

The  system  has  two  options  for  operation:  a)  Multi- 
line  operation,  b)  Single  line  operation.  The  output 
mirror  should  be  replaced  by  the  diffraction  grating  for 
single  line  operation,  while  two  mirrors  (one  partially 
transmitting)  are  used  for  multiline  operation. 

Optimized  power  output  from  lasing  of  hydrogen 
fluoride  was  reported  [Ref.  9]  as  6.5  Watts  with  the  mass 
flow  rates 

SFg  = 0.600  grams/sec. 

02  = 0.182  grams/sec. 

H2  = 0.019  grams/sec. 

He  = 0.053  grams/sec. 

Optimized  multiline  power  output  of  deuterium  fluoride 
lasing  was  also  reported  [Ref.  9]  as  2.5  Watts  with  the 
mass  flow  rates 
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SFg  » 0.680  grams/sec. 

0^  = 0.270  grams/sec. 

D2  = 0.043  grams/sec. 

He  = 0.063  grams/sec. 

In  the  case  of  single  line  operation,  the  wavelengths 
and  corresponding  relative  intensities  of  fourteen  individual 
DF  and  fifteen  individual  HF  lasing  transitions  are  shown 
in  figure  3 and  figure  4. 

The  block  diagram  of  the  main  elements  of  the  system 
is  also  shown  as  figure  5. 
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II.  OPTICAL  PROPAGATION  THROUGH  THE  ATMOSPHERE 

A.  MOLECULAR  SCATTERING  COEFFICIENT 

The  simple  Lorentz-type  molecular  model  presented  here 
is  shown  in  figure  1.  It  is  assumed  nonionized,  nonpolar, 
isotropic,  linear  and  lightly  damped. 

The  model  contains  molecular  mass  residing  in  the 
atomic  nuclei  and  carrying  a net  positive  charge.  Orbital 
electrons  shared  by  the  constituent  atoms  are  represented 
by  a concentric  shell,  providing  a balancing  negative 
charge.  Force  directed  toward  the  central  mass  is  char- 
acterized by  a spring  constant  k.  Let  m be  the  electron 
mass.  If  the  system  is  disturbed  by  an  electric  field  E 
in  the  z direction,  then  the  equation  of  motion  may  be 
written  for  a displaced  electron  as 

d^z 

- eE  - kz  = m 2-4  (2-1) 

z dt2 


where 


E = E sin  ait. 
z oz 


I 


Equation 


(2-1) 


can  also  be  written  as 


2 oz  • 

+ a)  z = sm  <i>t 

o m 


(2-2) 
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where  u)Q  = (k/m)  is  defined  as  the  resonance  frequency, 
corresponding  to  a molecular  transition. 

The  steady  state  solution  to  equation  (2-2)  is 


E sin  wt 


m(u)o  - u ) 


Now,  the  oscillator  dipole  moment  can  be  written 


p = -ez  = 


— =-  E sin  cot 

u,2  - u2  02 

o 


If  we  denote  the  maximum  value  of  the  dipole  moment 
as  pQ,  then  it  is 


-g2-ZP.  E 

2 2 boz 

U>  - (D 

o 


(2-3) 


The  oscillating  dipole  moment  has  maximum  value  pQ 
given  by  equation  (2-3)  and  produces  secondary  waves. 
Reradiation  from  the  oscillating  dipole  then  produces 
electric  and  magnetic  fields  at  every  point  in  the  coor- 
dinate system  shown  in  figure  2.  In  the  short  dipole 
approximation  (i.e.,  radius  of  scatterer  is  much  less  than 
wavelength  of  incident  wave)  and  far  field  (i.e.,  distance 
to  scatterer  is  much  larger  than  radius  of  scatterer)  the 
electric  and  magnetic  field  vectors  can  be  written  as 
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where 


E = — US:L---  sin  (kr-wt ) 9 


B = MS^"- 9 sin(kr-wt)  <fi 


v is  the  velocity  of  wave, 
y is  the  magnetic  permeability, 

0 , 0 and  r are  the  variables  in  spherical  coordinates 


The  poynting  vector  denoted  by  S describes  scattered  power 
per  unit  area,  and  is  given  here  as 


S = 


^(^£)2(£%^-)2sin2(kr-wt)^  (2-4) 


Scattered  power  per  unit  solid  angle  can  be  found  from 

equation  (2-4),  first  taking  the  time  average,  then  multi- 
2 

plying  by  r . This  gives  the  scattered  intensity  as 


- £^)2sin2(e>  = 


y , oi  e \2p2 

Tv  . ,2  2.  boz 

4iTm(oj  - w ) 
o 


(2-5) 


the  power  of  the  incident  wave  contained  in  unit  area  can 
be  given  as 


<S>incident 


\ v e E2 

Cm  O Z 


(2-6) 


where  e is  the  permittivity  of  the  medium. 
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The  angular  cross  section  of  a molecule,  denoted  by 
am(9),  may  be  defined  as  that  cross  section  of  incident 
wave  acted  on  by  the  molecule,  having  an  area  such  that 
power  flowing  across  it  is  equal  to  the  power  scattered 
by  the  molecule  per  unit  solid  angle  at  an  angle  9.  It 
is 


a (9)  = 

m 


a o) 
m 


^scat(8) 

<S>incident 


equation  (2-5) 
equation  (2-6) 


or 


' 4 7T  ' ' 


*■  M 

e x 2 

2 T~ 
m(w_  - u ) 


/Uco  , 2 . po  . 2 

1 4 IT  ' (E  ' 
oz 


(2-7) 


It  is  necessary  to  define  the  polarization  P in  order  to 
find  a suitable  form  for  the  second  term  in  the  right  hand 
side  of  equation  (2-7) . P is  defined  as  total  dipole 
moment  per  unit  volume  and  it  is  given  as 


p = Npo  = XeoEeff  = (K-l)€0Eaf£ 


for  a linear  and  isotropic  medium,  where 

N is  number  of  molecules  per  unit  volume, 

X is  electric  susceptibility, 

K is  dielectric  constant, 

is  the  effective  electric  field  which  is 

E + . 

oz  3e_ 
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Now  we  can  write  P in  terms  of  K and  E , as 

oz 


3 (K-l) e 

2 p 

K+2  oz  ' 


from  which 


oz 


3(K-l)e{ 

“nTkTTT 


The  index  of  refraction  is  given  here  as 


n = § « SOT  = /K 

v m 


since  Km  = 1 for  a nonmagnetic  medium. 
Now  equation  (2-7)  can  be  written  as 


...  2 3e  ( n — 1 ) 

°m!e>  = <TF>(-2T > 

N(n  +2) 


defined  by  the  equation  below,  is  called  the  total 
scattering  cross  section. 


or 


4tt 

J = / a (0 ) dft 

m ' m 


8tt3  / 3 (n2-l)  \ 2 


m 


2 4 v 2 
3N^X  (n  +2) 


(2-8) 
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where 

2itv 

“ = T- 

was  substituted. 

The  total  scattering  cross  section  given  by  equation 

(2-8)  is  defined  similarly  to  the  angular  scattering  cross 

section,  as  the  power  scattered  in  all  directions. 

For  air  the  refractive  index  is  very  close  to  one. 

2 

Therefore  n +2  is  nearly  equal  to  three,  and  equation  (2-8) 
may  be  written  as 


m 


8ttj  , 2 .,2 
(n  -1) 


3N2  A4 


(2-9) 


This  form  of  the  total  scattering  cross  section  (2-9) 
has  been  derived  for  a single  molecule  disturbed  by  incoming 
electromagnetic  radiation.  It  has  the  units  of  area.  For 
many  molecules  in  a unit  volume,  the  total  scattering 
coefficient  Sm  may  be  defined,  assuming  each  molecule  has 
the  same  total  scattering  cross  section, 


m 


= No- 


rn 


5V(n2-l) 


3NA 


(2-10) 


The  molecular  scattering  coefficient  (2-10)  has  the  units 
of  1/length. 

Since  atmospheric  gases  are  not  quite  isotropic,  a 
depolarization  factor  is  introduced  [Ref.  1]  and  inserted 
into  equation  (2-10) . The  molecular  scattering  coefficient 
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may  then  be  written  as 


6 


m 


8 j(n2-l) 2 


3NA 


+ 3Pd 
§ - 7P 


taking  into  account  the  pressure  and  temperature  dependence. 

Here  P_  and  T refer  to  standard  conditions,  1013  mb, 
o o 

273  K. 

Assuming  air  to  be  an  ideal  gas  and  using  the  ideal 
gas  law,  air  is  found  to  have  a mole  density  of  44.631  moles 
per  cubic  meter  at  a pressure  of  1013  millibars  and  tempera- 
ture of  273  degrees  Kelvin.  This  is  equivalent  to 
25 

2.688  * 10  molecules  per  cubic  meter. 

For  a depolarization  factor,  a value  of  0.035  is 
used,  following  Penndorf  [Ref.  1] . Although  the  index  of 

refraction  of  air  depends  on  the  wavelength,  this  variation 

2 2 -7 

is  small.  Therefore  (n  -1)  is  taken  as  3.0884  x 10 

which  corresponds  to  a wavelength  0.55  micrometers  [Ref.  2]. 

Using  these  parameter  values  given  above,  the  atmospheric 

molecular  scattering  coefficient  may  be  written  as 


B 


m 


10.07  x 10-20  , P , ,273, 
'1013' ' T ' 


units  of  1/km 


where 

A has  units  of  cm. 

P has  units  of  mbars. 

T has  units  of  degrees  Kelvin. 


McClatchey  and  d'Agati  [Ref.  4]  have  used  the  slightly 


different  form  of  the  molecular  scattering  coefficient, 


9.807  x 10“20  , P , ,273, 
~ ^4 . 0117  ^1013  1 T 1 


(2-11) 


which  was  selected  as  giving  a "best  fit"  to  the  experi- 
mental data  of  Penndorf  [Ref.  6]. 

B.  MOLECULAR  ABSORPTION  COEFFICIENT 

The  model  to  derive  the  molecular  absorption  coeffi- 
cient is  the  same  as  the  one  we  developed  in  the  previous 
section,  but  now  it  is  heavily  damped.  When  the  model  is 
disturbed  by  an  electric  field,  the  electron  moves  with  a 
damping  coefficient  and  restoring  force  strength  k.  The 
equation  of  motion  may  be  written  in  differential  form  such 
that 


m + my  + kr  = 
dt2  dt 


-eE  = -eE(r)e 


-loot 


The  steady  state  solution  to  this  equation  is 


r = 


-e/m 


2 2 

go  - go  - lyu) 
o 


The  polarization  due  to  N molecules  per  unit  volume  is 


P = -eNr  = 


e2N 


m 


1 

T 


CO  — GO  - lyco 


'eff 


(2-12) 
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for  the  index  of  refraction. 

A more  general  form  for  equation  (2-14)  is  written 
below  with  j referring  to  different  types  of  molecules 
with  oscillator  strength  f y 


Near  a resonance  frequency,  equation  (2-15)  can  be 
approximated  by 


e2N  . f 


[1  + 


2_1 


2e  m 
o 


. 2 2\ 


, 2 2.  2 , 2 2 
(oj1  j— a)  ) +y  w 


1 + c 


e2N  f . 


2e  m 
o 


-“Z- 


2 2 . 2 2 2 
UJlj-u  ) +Y  u> 


] i 


(2-16) 


using  the  Binomial  Expansion,  since  n is  very  close  to 
unity. 

A plane  wave  going  in  the  z direction  with  wave  vector 
K may  be  written  as 


E - I ei(K*-"t) 
o 


(2-17) 


The  relation  between  wave  vector  and  refractive  index 
is 


K = — n (2-18) 

c 

K is  complex,  since  n is  complex,  so  defining  K and  n 

K = 6 + ia  (2-19) 

n = n + ik  (2-20) 

where  S,  a,  n,  k are  the  so-called  propagation  constant, 
absorption  coefficient,  index  of  refraction,  and  extinction 
coefficient  respectively. 
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Using  equations  (2-16,18,19,20),  it  can  be  shown  that 


e2^ 


2 2 
a)  Y 


a = 2e  me,  Y ~2  2.2^  2 2 

° / (wlj-ai  ) +oj  Y 


(2-21) 


Substituting  K from  equation  (2-19)  into  equation  (2-17) 
one  can  get 


E = E e"az  ei(8z-wt) 
o 


The  electric  field  decreases  by  a factor  of  e-az, 

• “ 2 ft  2 

and  the  intensity  falls  off  by  a factor  of  e along  the 
path  length  z.  Since  a corresponds  to  absorption,  many 


^ ¥ 


where  A is  a constant,  Av  is  y/4ir  and  defined  as  half 
width  at  half  maximum. 

If  a is  normalized  in  the  range  from  minus  infinity  to 
plus  infinity,  the  expression  for  the  Lorentzian  line 
shape  becomes 


g(v) 


Av 


it  [ (vQ-v)  2 + Av^] 


describing  the  spectral  distribution  of  absorption. 

Finally,  the  molecular  absorption  coefficient  can  be 
written  as 


K (v) 
m 


Set 


2 T~  ' 

7T  [ (vQ-v)  + a ] 


provided  S = /Km(v)dv  for  single  absorbing  line,  where 


S is  line  intensity  per  absorbing  molecule  in 
the  units  of  cm~V  (molecules/cm^)  , 


a is  Av  half  width  at  half  maximum  in  the  units 
of  cm-1, 


-1 


v is  incoming  wave  frequency  in  the  units  of  cm  , 


; is  resonance  frequency  for  absorbing  line  in 
the  units  of  cm 


The  radiant  frequency  v is  unique:  this  gives  absorption 
as  a spectral  quantity.  Radiation  from  laser  sources  has 
very  narrow  line  width,  which  may  be  assumed  infinitely 
narrow  for  calculation  purposes.  The  resonance  frequency 
for  each  absorbing  line  is  independent  of  pressure  and 
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temperature.  The  line  width  is  assumed  to  be  a pressure 

broadened  Lorentzian,  with  its  pressure  and  temperature 

dependence  given  by  the  expression  [Ref.  7] 

T 

a = “(W  (~T)1/2  (2-23) 

o 


where 


P_  = 1 atm.  and 

o 

Tg  = 296  degrees  Kelvin. 

The  total  intensity  of  any  line  is  pressure  independent 
and  its  temperature  dependence  is  given  by  the  expression 
[Ref.  7] 


QV(T  ) Q (T  ) T-T 

S(T)  = S ( Tg ) q^(tT~  * Q~TtT  * exp  [ 1 • 4 3 9 E " ( ) ] (2-24) 

where  Qv  and  Qr  are  vibrational  and  rotational  partition 
functions  respectively,  E"  is  the  lower  state  energy  of 
the  transition  corresponding  to  that  line,  in  units  of  cm  1. 

Rotational  partition  functions  are  also  temperature 
dependent,  of  the  form 


w 

Qr(T) 


j 


The  corresponding  values  of  Qv  and  j are  tabulated  for  each 
molecule  in  the  air  [Ref.  7]  and  presented  in  table  2. 
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Finally  the  total  molecular  absorption  coefficient  due 
to  water  vapor,  ozone,  carbon  dioxide,  methane,  carbon 
monoxide,  oxygen,  and  nitrous  oxide  in  the  air  is  given 
by 


K. 


m 


l l 

j i 


S. .a. .m. 


ir[  (v-vij  ) 2 + ajj] 


(2-25) 


Summation  over  j corresponds  to  different  types  of  mole- 
cules, and  summation  over  i corresponds  to  different  lines 
for  a given  molecule. 

Values  of  vj_j/  ^ij  (PQ,Tg)  , (Tg)  have  been  compiled 
for  100,000  lines  between  1 ym  and  the  far  infrared  by 
McClatchey  and  others  at  Air  Force  Cambridge  Research 
Laboratories  and  are  available  on  magnetic  tape  [Ref.  7] . 

Also  a new  version  of  the  compilation  including  over 
139,000  lines  between  0.2  ym  and  30  ym  was  recently  reported 
by  Air  Force  Geophysics  Laboratory  [Ref.  3] . 

Calculation  for  Km  at  any  pressure  and  temperature  may 
be  made  by  using  equations  (2-23,24,25)  with  the  data  from 
this  tape. 

nr  in  equation  (2-25)  represents  the  number  of  j type 

molecules  per  unit  volume.  The  resultant  K is  in  the  units 

m 

of  1/length  and  generally  1/kilometer  is  used. 

It  should  be  noted  that  the  Lorentzian  line  shape  must 
be  modified  for  pressures  less  than  100  millibars;  at  pressures 
between  10  and  100  millibars  doppler  broadening  becomes 


1 
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dominant  [Ref.  7],  while  at  pressures  less  than  10  milli- 
bars, the  Voigt  profile  should  be  used.  The  Lorentzian 
line  shape  is  good  enough  to  be  used  at  low  altitudes,  but 
line  wings  should  be  truncated  at  20  cm  1 from  the  line 
centers,  since  uncertainties  begin  to  arise  beyond  that 
range. 

A computer  program  calculating  the  molecular  absorption 
coefficient  is  prepared  and  presented  as  Appendix  A.  This 
program  also  includes  the  water  vapor  continuum  absorption 
calculation  which  will  be  explained  in  the  next  section. 

The  algorithm  for  the  program  is  also  discussed  in  Section 
III. 


C.  CONTINUUM  ABSORPTION  COEFFICIENT 

Continuous  regions  of  absorption  in  the  atmospheric 
windows  occur  near  2500  cm  1 (4  ym)  and  from  700  cm  1 to 
1250  cm"1  (8-14  ym)  . 

In  other  spectral  regions,  the  contribution  from  nearby 
absorption  lines  is  much  greater  than  that  from  the  continuum 
absorption,  so  that  for  practical  purposes  the  continuum 
effect  can  be  neglected,  although  it  may  be  greater  than 
in  the  windows . 

Extreme  wings  of  strong  collision  broadened  absorption 
lines  centered  more  than  20  cm  1 away  or  pressure  induced 
absorption  resulting  from  transitions  that  are  forbidden  for 
unperturbed  molecules  may  cause  continuous  absorption.  In 


the  8-14  pm  region,  existence  of  the  water  dimer  (H20:H20) 
may  also  cause  continuous  absorption. 

1.  Water  Vapor  Continuum 

The  absorption  coefficient  due  to  the  water  vapor 
continuum  can  be  described  empirically  by  the  expression 
[Ref.  4] 

K = [Cj,  ( v , T)  P + C.T  (v  ,T)  P J W (2-26) 

mcont.  S s N N 

where  Cg  is  a self  broadening  coefficient  due  to  self 
broadening  from  collisions  of  water  molecules  with  other 
water  molecules,  CN  is  a foreign  gas  broadening  coefficient 
due  to  nitrogen  broadening  from  collisions  of  water  mole- 
cules with  air  molecules  which  are  mostly  nitrogen. 

Cg  and  CN  are  both  temperature  and  frequency  dependent, 
a.  8-14  Mm  region  [Ref.  4] 

CN  is  given  by  the  expression 

Cn(v,T)  = 0.002  x Cg ( v , 296 ) . 

is  given  by  the  expression 

Cg ( v , T)  = Cg ( v , 296 ) x exp [1800  (^  - ] 

where  Cg ( v ,296)  is  given  as 
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CQ(v,296)  - 4.18  + 5578  * exp(-7.87  x lo"J  v) 


where  Cg  is  in  the  units  of  (precipitable-cm) -1  atm  1 or 
(pr.  cm)  1 atm  1 which  is  (grams/cm2) -1  atm-1, 
b.  3. 5-4. 2 ym  region  [Ref.  4] 

C-  and  CM  are  given  by  the  expressions  below; 


Cn(v,T)  = 0 . 12Cg ( v , T) 


Cg(v,T) 


Cg  (v  , 296 ) x exp  [1350  (^  - ^)  ] 


where  Cg  has  the  units  of  (pr.  cm)  1 atm  1.  The  values  of 
Cg(v,296)  at  different  frequencies  are  reproduced  [Ref.  4] 
and  tabulated  in  table  3.  W in  the  equation  (2-26)  has 
the  units  of  (pr.  cm)  per  kilometer  which  is  the  same  as 
grams  per  square  centimer  per  kilometer.  Then  the  water 


vapor  continuum  absorption  coefficient 
obtained  in  (kilometers)  1. 


can  be 


cont . 


2 . Nitrogen  Continuum 

Although  the  nitrogen  continuum  is  not  included  in 
this  thesis,  the  magnitude  of  this  effect  is  of  the  order 
of  0.001  km  1 at  a wavelength  of  3.8  ym.  The  maximum 
effect  which  is  of  the  order  of  0.1  km  1 occurs  at  a wave- 
length of  4.3  ym.  This  effect  may  be  neglected  outside 
the  3. 7-4. 7 ym  region  [Ref.  4]. 
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where  r is  the  radius  of  the  scatterer  or  absorber,  and  X 
is  the  wavelength  of  the  incident  radiation,  a is  the 
parameter  on  which  aerosol  scattering  and  absorption  depend 
most  strongly,  rather  than  on  radius  of  particle  alone  or 
on  the  wavelength  of  incident  radiation  alone.  The  refrac- 
tive indices  used  in  calculations  are  the  relative  refrac- 
tive indices  of  the  aerosol  particles  with  respect  to  the 
refractive  index  of  the  medium  surrounding  the  particles. 
Since  air  has  refractive  index  very  close  to  unity,  the 
relative  refractive  index  is  assumed  equal  to  the  absolute 
refractive  index  of  the  aerosol  particle,  and  is  generally 
represented  by  a complex  number,  where  the  imaginary  part 
takes  care  of  possible  absorption  by  the  particle.  The 
refractive  indices  of  different  kinds  of  particles  used  in 
calculations  are  reproduced  and  presented  in  table  10 
[Ref.  12],  taking  into  account  the  wavelength  dependence. 

Total  light  scattered  or  scattered  and  absorbed  by  a 
single  particle  is  given  by  its  total  scattering  or 
extinction  cross  section  which  is  related  to  the  geometric 
cross  section  by  a parameter  generally  called  the  effi- 
ciency factor.  According  to  the  Mie  theory,  the  scattering 
and  extinction  efficiency  factors  are  given  by  the  equations 


l (2n+l) [|an|2+|bn|2] 
n=l 


(2-28) 


Thus  the  efficiency  factor  is  defined  as  the  ratio  of 
apparent  cross  section  to  geometric  cross  section  of  the 
particle.  a , a in  the  above  are  the  apparent  scattering 

S C 6X 


and  extinction  total  cross  sections  respectively.  Re 

represents  the  real  part  of  the  argument.  The  coefficients 

an,  bn  are  the  complex  functions  which  represent  the 

th 

amplitude  of  n electric  partial  wave  and  n magnetic 
partial  wave  respectively.  These  coefficients  may  be 
found  from  Ricatti-Bessel  functions  for  a given  complex 
refractive  index  and  size  parameter  a. 

The  scattering  efficiency  factor  approaches  the  value 
of  two  in  the  short  wavelength  approximation,  which  implies 
that  the  scattering  cross  section  is  almost  twice  as  large 
as  the  geometric  cross  section  of  the  particle. 

Since  the  infinite  series  represented  by  equations 
(2-28)  and  (2-29)  are  slowly  converging,  calculations  are 
generally  computer  time  consuming.  For  particles  with  real 
part  of  refractive  index  between  1 and  1.5  and  imaginary 
part  of  refractive  index  between  0 and  0.25,  the  Van  de 
Hulst  approximation  may  be  used  with  some  correction  fac- 
tors. Within  the  approximation  range,  it  was  reported  by 
Deirmendjian  [Ref.  13]  that  the  efficiency  factor  could 
be  predicted  within  the  precision  range  of  ± 0.05  Q . 

The  method  of  calculation  using  the  Van  de  Hulst 
approximation  is  presented  in  Appendix  C. 

The  efficiency  factor  of  a single  particle  then  may 
be  extended  for  use  in  calculation  of  the  aerosol  extinction 

I 


j 


iT,rr- =• 


coefficient  due  to  many  particles  for  a given  size 
distribution  such  that 


U 


a 


-3  / a2  Qsc  n (a)  da 

" *1 


a. 


-3  / a2  Qav  n(a)  da 


•ex 


a. 


(2-30) 


(2-31) 


where  K is  the  propagation  constant  and  n(a)  is  the  size 
distribution  defined  by  equations  (1-6)  and  (2-27) . 

Once  the  aerosol  scattering  and  extinction  coefficients 
are  determined,  then  the  aerosol  absorption  coefficient  may 
be  found  from  K = u - 8 , if  needed. 

dad 

Size  distributions  used  here  for  calculation  of  aerosol 
extinction  coefficient  are  a)  equation  (l-6a)  which  is 
independent  of  wind  velocity  and  relative  humidity  having 
75%  sea  spray  and  25%  rural  aerosol  mixture  which  is  com- 
posed of  30%  dustlike  and  70%  water  soluble  substance, 
b)  equation  (l-6b)  which  adds  the  effects  of  wind  velocity 
and  relative  humidity. 

If  the  total  number  of  particles  is  not  known,  then  the 
size  distribution  n(a)  or  n(r)  is  generally  normalized  to 
one  particle  per  cubic  centimeter.  Whatever  value  the 
size  distribution  is  normalized  to  we  can  put  a normaliza- 
tion constant  into  equations  (2-30)  and  (2-31) . Thus  we 
may  write 
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6 

a 


-t  / 


a 


K' 


sc 


n (a) 


a. 


da 


where  C is  the  normalisation  constant  and  may  be  found 
from  the  relation  [Ref.  19] 


R ~ B (0 ! 55  pm) 

CL 

Here  R is  the  visibility  in  kilometers  providing  8 (0.55) 

Ci 

in  inverse  kilometers. 

Therefore  knowing  the  range  R for  visibility,  the  aero- 
sol scattering  coefficient  can  be  found  from  the  equation 
above.  The  aerosol  scattering  coefficient  at  0.55  micro- 
meters calculated  by  equation  (2-30)  then  allows  determination 
of  the  value  of  the  normalization  constant  C.  Once  C is 
determined  at  0.55  micrometers,  it  is  the  same  for  the 
whole  range  of  wavelengths . 

Since  aerosol  absorption  by  particles  is  much  less  than 
scattering,  the  visibility  range  may  also  be  written  as 

r = 3.9 

U (0.55  ym) 

Si 

The  curves  of  extinction  coefficient  versus  wavelength 
for  varied  wind  speed  and  relative  humidity  shown  as  figures 
10  and  11  are  normalized  so  that  the  total  aerosol  extinction 
coefficient  at  0.55  um  is  0.17  km  ^ . This  corresponds  to 
a visibility  of  23  kilometers. 
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III.  ALGORITHM  FOR  CALCULATION  OF 

MOLECULAR  ABSORPTION  COEFFICIENT 


The  computer  program  presented  as  Appendix  A is 
intended  to  calculate  the  molecular  absorption  coefficient 
including  the  water  vapor  continuum  effect.  As  described 
in  Section  II. B,  the  program  reads  the  line  parameters 
(Vij,  cu  j , Sij)  from  a magnetic  tape.  However,  other 
inputs  of  the  program  are  the  amounts  of  the  air  molecules 
which  have  to  be  specified  in  the  units  of  molecules  per 
square  centimeter  per  kilometer  path,  and  they  should  be 
supplied  by  the  user. 

The  amounts  of  different  types  of  air  molecules  excluding 
the  ozone  and  water  vapor  are  calculated  from  the  ideal  gas 
law. 


P = p R T 

where  P is  air  pressure  in  Newtons  per  square  meter,  T 
is  air  temperature  in  degrees  Kelvin,  R is  the  gas  con- 
stant for  air  in  Joules  per  kilogram  per  degree  Kelvin 
and  has  the  value  of  287.06  Joules  per  kilogram  per  degree 
Kelvin.  For  given  values  of  pressure  and  temperature,  the 
density  of  air  can  be  found  in  the  units  of  kilograms  per 
cubic  meter.  Then  the  density  of  air  is  used  to  find  the 
amounts  of  the  uniformly  mixed  gases  (02,  N2,  CC>2,  CO, 

N20,  CH4)  such  that 


t 


I I 
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p.  = x.  10  p . 

i i air 


where  is  the  concentration  in  PPM  by  volume  taken  from 
table  4,  and  p^  is  the  density  or  the  mass  per  unit  volume 
of  the  corresponding  constituent.  After  finding  the  amount 
of  a given  type  of  molecule  in  air  in  the  units  of  mass 
per  unit  volume,  it  is  necessary  to  convert  it  to  the 
units  of  molecules  per  unit  volume  to  be  able  to  use  it 
in  the  computer  program.  This  can  be  done  through  the 
relation 


p . N 

w = x n 
i -3  ' 

Mi  10  J 


where  p^  has  the  units  of  kilograms  per  cubic  meter, 

+23 

is  Avagadro ' s number  which  is  6.023  * 10  molecules  per 
mole  and  is  the  molecular  weight  in  grams  per  mole  of 
the  corresponding  constituent.  Now  the  amount  of  the  ith 
constituent  of  air  is  given  by  the  formula  above  in  units 
of  molecules  per  cubic  meter.  Final  conversion  may  be 
made  into  units  of  molecules  per  square  centimeter  per 
kilometer  which  may  be  used  directly  as  input  to  the  com- 
puter program  presented  in  Appendix  A.  may  also  be 

expressed  in  terms  of  P and  T as 


19  P Xi 

Wi  = 2.0982  x ioiy  £ jji 

i 


(3-1) 


where  P is  air  pressure  in  millibars,  T is  temperature  in 

degrees  Kelvin,  X^  and  are  as  defined  before.  Now 

is  the  amount  in  molecules  per  square  centimeter  per 

kilometer  for  the  ith  type  of  air  molecules. 

Using  the  equation  (3-1) , the  amount  of  carbon  dioxide 

20 

which  has  X^  = 330  and  = 44  is  5.5943  * 10  molecules 
per  square  centimeter  per  kilometer  in  the  air  where 


P = 1017.5  millibars  and  T = 286.22  degrees  Kelvin. 

The  densities  of  non  uniformly  mixed  gases  are  more 
complicated  to  determine.  For  ozone  at  sea  level,  it  is 

_ 5 

taken  as  6 * 10  grams  per  cubic  meters  or  7.5288  x 10 
molecules  per  square  centimeter  per  kilometer  at  any  given 
temperature  and  pressure. 

The  amount  of  water  vapor  and  its  partial  pressure  can 
be  found  from  the  dew  point  temperature  or  relative  humidity 
in  the  air  [Ref.  15].  Different  dew  point  temperatures 
and  corresponding  density  and  partial  pressure  of  water 
vapor  are  tabulated  in  table  11. 

The  water  vapor  density  determined  by  interpolating 
from  table  11  may  be  converted  into  number  density  (per 
square  centimeter  per  kilometer  path)  as  previously 
described. 


IV.  RESULTS 


Although  several  band  models  were  introduced  by  many 
workers  in  order  to  fit  the  molecular  absorption  lines  in 
the  atmosphere,  line  by  line  calculations  are  made  here, 
making  use  of  the  experimentally  found  and  compiled  data. 

The  molecular  absorption  coefficient  is  the  most  important 
attenuation  coefficient  to  be  calculated,  since  mainly  it 
is  the  one  which  determines  the  atmospheric  windows  which 
are  good  enough  for  propagation  of  light.  Also  it  depends 
strongly  on  the  frequency  of  propagation.  For  any  fre- 
quency of  laser  radiation,  the  molecular  absorption  coeffi- 
cient can  not  be  computed  by  interpolating  between  two  close, 
known  values.  It  can  be  computed  either  by  the  line  by 
line  method  (infinite  resolution)  or  by  interpolation  from 
calculations  with  resolution  less  than  0.01  cm-1.  Therefore 
the  reported  values  of  molecular  absorption  coefficients 
may  be  used  directly.  If  any  other  frequency  or  condition 
is  needed,  the  computer  program  presented  as  Appendix  A may 
be  used.  Atmospheric  windows  in  the  3-5  micrometer  and  8-14 
micrometer  regions  are  mainly  available  for  atmospheric 
transmission.  HF/DF  laser  wavelengths  fall  in  the  near 
and  middle  infrared  region  of  the  spectrum.  In  particular 
the  P2 ^ transition  of  the  DF  laser  with  wavelength  3.8007 
micrometers  shows  relatively  high  transmittance  falling 
in  the  3-5  micrometer  region  of  the  spectrum.  The  molecular 


absorption  coefficient  for  this  particular  wavelength  is 
0.0104  km  1 in  the  midlatitude  winter  model  where  the 
molecular  scattering  coefficient  is  essentially  zero. 

This  shows  99  percent  transmittance  along  a one  kilometer 
path. 

Calculations  intended  for  use  in  the  propagation  studies 
in  the  Monterey  Bay  area  make  use  of  the  known  range  of 
conditions  in  Monterey,  taking  into  account  the  yearly  and 
monthly  averages  of  air  pressure,  air  temperature  and  dew 
point  temperature  in  the  Monterey  Bay  area.  These  yearly 
and  monthly  averages  are  taken  from  [Ref.  16]  and  [Ref.  17] ; 
the  quantities  used  in  calculations  are  tabulated  in  table 
12.  In  this  table  the  mean  temperature  is  used  which  is 
the  mean  value  of  maximum  and  minimum  average  temperatures 
corresponding  to  day/night  condition. 

Also  the  amounts  of  air  constituents  used  in  the  calcu- 
lations are  reported  in  table  17. 

Table  9 shows  the  calculated  molecular  absorption  for 
some  selected  wavelengths  for  the  midlatitude  winter  model, 
using  the  computer  program  presented  in  Appendix  A.  Table 
13  contains  the  molecular  absorption  coefficients  for  yearly 
averages  in  the  Monterey  Bay  area. 

Tables  14,  15,  16  show  the  molecular  scattering  and 
absorption  at  1.06,  3.8007  and  10.591033  micrometers 
respectively  for  the  monthly  averages  in  the  Bay  area. 

Since  the  molecular  scattering  is  negligible,  it  is  assumed 
to  be  zero  at  3.8007  and  10.591033  micrometers. 
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The  aerosol  extinction  coefficients  are  calculated 
for  the  rural  and  maritime  models  using  the  size  distri- 
bution given  by  equation  (l-6a)  and  presented  in  figure  8 
and  figure  9 respectively.  The  maritime  model  with  75% 
sea  spray  is  normalized  so  that  the  total  number  of  parti- 
cles is  about  44  per  cubic  meter.  These  results  agree  with 
the  previously  reported  values  [Ref.  12]  within  1.8  percent. 

The  relative  humidity  and  wind  velocity  dependent  size 
distribution  is  used  in  figure  10  and  figure  11,  normalizing 
the  curves  to  a visibility  of  23  kilometers.  In  both 
figures  71%  sea  spray  and  29%  rural  aerosols  are  assumed. 
Although  at  0.55  micrometers,  the  rural  aerosol  refractive 
indices  are  outside  the  approximation  range,  the  results 
seem  to  be  valid  since  the  percentage  of  rural  aerosols 
is  low. 

As  seen  in  figure  10,  increasing  wind  speed  increases 
the  attenuation  in  the  longer  wavelength  region.  Similarly 
in  figure  11,  increasing  relative  humidity  increases  the 
attenuation  in  the  longer  wavelength  region. 

Figures  12  through  19  show  the  transmittance  along  a 
path  length  of  10  kilometers , due  to  the  molecular  absorption 
only. 

Figures  12,  14,  16,  18  are  produced,  using  a spectral 
resolution  of  0.01  cm  1 for  use  with  laser  sources. 

A degraded  resolution  is  used  in  figures  13,  15,  17,  19 
by  averaging  the  transmittance  at  each  point  over  a ?* ectral 
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bandwidth  of  0.5  cm  Therefore  the  figures  13,  15,  17, 

19  are  for  use  with  the  broadband  sources  only. 

Each  one  of  the  figures  12,  13,  16,  17,  18,  19  covers 
a spectral  interval  of  about  25  cm  1 near  the  wavelength  of 
3.8  micrometers  (i.e.,  2631. 09. cm  1) . The  figures  14  and 
15  are  for  the  region  about  the  wavelength  of  10.6  micrometers 
(i.e.  , 943.936  cm"1) . 

Figures  12  through  15  show  the  transmittances  for 
yearly  averages  in  the  Monterey  Bay  area. 

Figures  16  through  19  are  produced  for  the  different 
partial  pressures  of  water  vapor,  having  the  same  air 
temperature  of  290  degrees  Kelvin  and  the  air  pressure  of 
1018  millibars.  The  partial  pressure  of  water  vapor  in 
figures  16  and  17  is  7.48  mbars,  and  it  is  15.48  mbars  in 
figures  18  and  19. 

Using  the  high  resolution  of  0.01  cm  1 and  a degraded 
resolution,  a sensitivity  analysis  due  to  the  air  pressure, 
temperature  and  partial  pressure  of  water  vapor  has  been 
carried  out  near  the  wavelength  of  3.8  micrometers. 

It  has  been  observed  that  the  change  in  the  air 
pressure  and/or  temperature  which  could  affect  the  line 
widths,  line  intensities  and  amounts  of  the  absorbing  air 
molecules,  does  not  contribute  to  the  transmittance  signi- 
ficantly, at  the  region  near  3.8  micrometers.  The  reason 
could  be  the  small  number  of  absorbing  lines  which  are 
mostly  due  to  water  vapor.  Comparing  figures  12,  16,  18  or 
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appendix  c 


CALCULATION  OF  EXTINCTION  EFFICIENCY  BY  USING 
THE  VAN  DE  HULST  APPROXIMATION 

This  approximation  was  first  developed  by  Van  de  Hulst 
[Ref.  13]  following  the  Huygens  principle.  It  is  good  for 
spherical  particles  of  arbitrary  size  but  |m|  •+  1. 

The  refractive  index  of  the  particle  is  denoted  by  m 
which  has  real  and  imaginary  parts. 

m = n - in . 

r i 

Defining  the  normalized  size  parameter  p and  absorption 
parameter  g such  that 


P 


g 


2a (nr  - 1) 


the  extinction  efficiency  may  be  written  as 


Q'ex  = 2 - - exp[~ptan  g]  sin  (p-g) 

+ 4 (£Os_£) 2 {cos  2g  - exp[-ptan  g]  cos(p-2g)} 

This  uncorrected  extinction  efficiency  Q^x  overestimates 
the  cross  section  for  small  p and  underestimates  it  in 
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varying  degrees  as  p approaches  and  surpasses  the  value 
4.08,  corresponding  to  the  first  maximum  in  Q'  . 

6X 

A correction  to  the  extinction  efficiency  is  made  in 
the  form 


‘ex 


(1  + D) 


where  D,  the  correction  factor,  is  given  by  the  expressions 


(n  -1) 


0.2p-n  +1 


D = 


= r:6«~n.lt(g)+11  * ir-rifigi  when  p-5(nr'l) 


D = i7ir4'  lf(g)tllp  when  5<nr-l>  1 P < l.kan 


(n  -1) 


D = 


" r “•  rj=  , , , , , , 4.08  , 4. OB 

2n,(l+3tan  g)  [£(g)+1]  when  I+3tirTg  - p - I+tan 


4.08 


D = 


2.04(nr-l) 
nrpf (g) 


■[f(g)+l]  when  p > 


1+tan  g 


where 


f (g) 


1 + 4 


tan  g + 3 (tan  g) 


2 


Within  the  limits  of 


1 < nr  <_  i.S 
0 < n . < 0.25 

— i — 
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Qgx  is  accurate  within  an  error  of  about  ± 0.05  Qgx 

for  a wide  range  of  sizes  and  types  of  index  of  refraction. 


TABLE  1.  MODIFICATION  FACTOR  FOR  C02  LINES 


V - V 

o 

X 

1/CM 

o 

• 

o 

1.00 

0.5 

1.00 

0.6 

0.96 

0.7 

0.89 

0.8 

0.32 

0.9 

0.77 

1.0 

0.70 

1.2 

0.60 

1.5 

0.50 

2.0 

0.  41 

2.5 

0.34 

o 

• 

CO 

0.31 

5.  C 

0.29 

00 

• 

o 

0.23 

10.0 

0.19 

15.0 

0.00 
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H20  1.5  1.0000  1.0000  1.0000  1.0000  1.0000  1.000  1.0000 
C02  1.0  1.0095  1.0192  1.0327  1.0502  1.0719  1.093  1.1269 
03  1.5  1.0040  1.0070  1.0130  1.0220  1.0330  1.046  1.0660 
N20  1.0  1.0170  1.0300  1.0480  1.0720  1.1000  1.127  1.1700 
CO  1.0  1.0000  1.0000  1.0000  1.0000  1.0000  1.000  1.0000 
CH4  1.5  1.0000  1.0000  1.0010  1.0020  1.0040  1.007  1.0110 
□2  l.c  1.0000  1.0000  1.0000  1.0000  1.0000  1.000  1.0010 


TABLE  3.  SELF  BROADENING  ABSORPTION  COEFFICIENTS  FOR 
WATER  VAPOR  IN  THE  3. 5-4. 2 MICROMETER  REGION 


V 

Cs ( v, 296 ) 

1/CM 

1/(PR.CM)ATM 

2350 

0.  230 

24-00 

0.187 

2450 

0.147 

2500 

0.117 

- 

2550 

0.097 

2600 

0.087 

2650 

0.100 

[ i 

2700 

0.120 

2750 

0.147 

2800 

0.174 

28  50 

0.200 

2900 

0.240 

2950 

0.280 

3000 

0.330 
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TABLE  4.  CONCENTRATIONS  OF  UNIFORMLY  MIXED  GASES 
IN  THE  AIR 


CONSTITUENT 


CONCENTRATION 
PPM  BY  VOLUME 


3 .300E+2 
2.300E-1 
7.500E-2 
1 •600E+0 
2.095E+5 
7.803E+5 


TABLE  5.  MODELS  OF  ATMOSPHERE  AT  SEA  LEVEL 


MODEL 


PRES. 

MB. 


TEMP.  DENSITY  H20 
DEG.  KEL.  GR/M3  GR/M3 


03 

GR/M3 


TROPICAL 


1013  300.0  1167  19.0  5.6E-5 


MIDLATITUDE 

SUMMER 

1013 

294.0 

1191 

14.0 

6.0E-5 

MIDLAT  ITUOE 

WINTER 

1013 

272.2 

1301 

3.5 

6.0E-5 

SUBARCTIC 

SUMMER 

1010 

287.0 

1220 

9.1 

4.9E-5 

SUBARCTIC 

1013 

257.1 

1372 

1.2 

4.1E-5 

WINTER 


U.S.  STANDARD  1013  283.1  1225  5.9  5.4E-5 

ATM.  ,1962 


IS  ■ ■ —I. — ; ~ 


TABLE  6.  SIZE  DISTRIBUTION  PARAMETERS  OF  AEROSOL  MODELS 


MODEL 


RURAL 


URBAN 


0.  9999975  0.005m  0.475  2.5E-6  0.5U  0.475 
0.9999975  0.005U  0.475  2.5E-6  0 . 5M  0.475 


MARITIME 

(CONTINENTAL) 

MARITIME 
(SEA  SPRAY) 


0.  005M  0.475 


0.3U  0.4 


N IS  NORMALIZED  TO  1 PARTICLE/CM3 


TABLE  7.  HF  LASER  FREQUENCIES  [10] 


VI 8RATIONAL 
TRANSITION  ID 


1-0 


2-1 


3-2  * 


4-3 


ROTATIONAL 
TRANSITION  ID 

P6 

P7 

P 8 

P9 

P10 

Pll 

P12 

P13 

PI  A 

P 15 

P2 

P 3 

PA 

P5 

P6 

P7 

PS 

P9 

P10 

Pll 

P12 

P13 

P1A 

P15 

P2 

P3 

PA 

P5 

P 6 

P 7 

P8 

P5 

P6 


EMISSION 

FREQUENCY  (1/CM) 

3693.50 
36AA.16 

3593.80 

3542.20 

3489.59 
3436.12 

3381.50 

3326.21 
3269.90 

3212.80 

3708.86 
3666.38 
3622.71 
3577.47 
3531.31 

3483.63 

3425.17 
3335.34 
3334.55 

3282.86 

3230.18 

3176.60 
3122.14 

3067.22 

3544.51 

3503.80 
3461.54 
3418.16 
3373.46 
3327.73 

3280.64 
3262.53 
3219.50 
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TABLE  7.  HF  LASER  FREQUENCIES  [10]  (Continued) 


5-4 


6-5 


P7 

3175.34 

P 8 

3130.09 

P9 

3083.83 

P4 

3150.67 

P 5 

3110.34 

P6 

3068.63 

P7 

3C26.21 

P8 

2982.51 

P9 

2937,79 

P5 

2961.68 

P 6 

2921.74 

P7 

2aeo.7o 

P8 

2838.59 
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If 


1 ' 1 

"*  ' 

TABLE  8.  DF  LASER  FREQUENCIES  [4] 

VI  6RAT  ION  At  ROTATIONAL 

EMISSION 

TRANSITION  ID  TRANSITION  ID 

FREQUENCY  (1/CM) 

r 

1-0  PI 

2 384.  934 

P2 

2862.652 

, • 

P3 

2839.779 

1 

P4 

2816.362 

P5 

2792.437 

P6 

2767.914 

P7 

2743.028 

P8 

2717.  536 

P 9 

2691.409 

P10 

2665.200 

Pll 

2633.396 

PX1 

2611.125 

P10 

2584.910 

P14 

2557.090 

P15 

2527.060 

* 

P 16 

2498.020 

2-1  P3 

2750.050 

P4 

2727.330 

' 

I ^ j 

P5 

27C3.980 

P6 

2680.280 

P 7 

2655.  970 

P 8 

2631.090 

P9 

2605.870 

P10 

2580.160 

i •’ 

Pll 

2553.970 

PH 

2527.470 

P13 

2500.320 

P16 

2417.270 

3-2  P3 

2662.170 

P4 

2640.040 

P5 

2617.410 

P6 

2594.230 

l| 

P7 

2570.510 

J 

TABLE  8. 

DF  LASER 

FREQUENCIES  [4]  (Continued) 

P8 

2546.370 

P9 

2521.810 

P10 

2496.610 

Pll 

2471.340 

P12 

2445.290 

P13 

2419.020 

P 14 

2392.460 

4-3 

P5 

2532.500 

P6 

2509.860 

P7 

2486.  830 

P8 

2463.250 

P9 

2439.290 

P10 

2414.890 

5-4 

P7 

2404.630 

7-6 

P8 

2222.680 

PIO 

2177.990 

PU 

2155.030 

P 12 

2131.680 

8-7 

P7 

2165.930 

P8 

2144.300 

P 9 

2123.240 

PIO 

2101.270 

Pll 

2056.140 

PIO 

2033.010 

9-8 

P6 

2108.480 

P7 

2038.340 

P 8 

2067.760 

PIO 

2025.360 

Pll 

2003.560 

P12 

1981.330 

77 


TABLE  9.  MIDLATITUDE  WINTER  MODEL  MOLECULAR  ABSORPTION 
COEFFICIENTS  FOR  INFINITE  RESOLUTION 


f 


FREQUENCY 
(1 /CM  I 


MOLECULAR  ABSORPTION 
COEFFICIENT  (I/KM) 


9 24*9  75 

4.27E-1 

942  .384 

6.17E-2 

944.195 

6.61E-2 

945.981 

6.66E-2 

9 47.743 

7.  22E-2 

959.393 

3.45E-2 

973.289 

7.49E-2 

974.623 

6.96E-2 

975.931 

1.32E-1 

977.215 

7.28E-2 

1046.854 

8.27E-2 

1048.661 

8.97E-2 

1050.441 

9.81E-2 

1075.988 

8.905-2 

1077.303 

8.21E-2 

1078.591 

7.37E-2 

2177.990 

3.89E-2 

2414.890 

1.16E-2 

2463.250 

1.26E-2 

2611.125 

7.19E-2 

2617.410 

7.19E-3 

2631.090 

1.04E-2 

2638.396 

7.90E-2 

2640.040 

1.67E-2 

2 655.97  0 

2.24E-2 

2662.170 

1.24E-2 

2665.200 

1 .64E-2 

2703.980 

9.55E-3 

2717.536 

3.3  7E-2 

2727.380 

1.38E-2 

2743.028 

1.16E-2 

2750.050 

1.61E-2 

TABLE  10.  REFRACTIVE  INDICES  OF  AEROSOL  PARTICLES  [Ref.  12] 


! 


WAVELENGTH 

WATER  SOLUBLE 

DUST 

LIKE 

SOCT  LIKE  OCEANIC 

(pm) 

nr 

ni 

nr 

ni 

nr 

ni 

n 

r 

n . 

l 

0.2 

1.53  0 

0.070 

1.53 

0.070 

1.  50 

0.35 

1.429 

2.87E-5 

0.25 

1.530 

0.  03  0 

1.  53 

0.030 

1.62 

0.45 

1.404 

1.45E-6 

0.3 

1. 53  0 

0.008 

1.53 

0.008 

1 .74 

0.47 

1 .395 

5.83E-7 

0.3371 

1.53  0 

0.005 

1.53 

0.008 

1 .75 

C.47 

1.392 

1.20E-7 

0.4 

1.530 

0.005 

1.  53 

0.008 

1.75 

C.  46 

1.3  85 

9.90E-9 

0.  4 6 8 

1.  530 

0.005 

1 . 53 

0.003 

1.75 

0.45 

1.382 

6.41E-9 

0.  5145 

1.530 

0.005 

1 .53 

0.008 

1.75 

0 .45 

1.381 

3.70E-9 

0.55 

1.53  0 

0 .006 

1 .53 

0.003 

1.75 

0.44 

1.  381 

4.26E-S 

0.6323 

1.53  0 

0.006 

1.53 

0.008 

i.  75 

0.43 

1.377 

1 .52E-3 

0.  6 943 

1.  530 

0.007 

1.53 

0.  008 

1.75 

0.43 

1 .376 

5 .04E-8 

0.36 

1 .520 

0.012 

1.52 

0.008 

1.75 

0 .43 

1.372 

1 .09E-6 

1.0  6 

1.520 

0.017 

1.52 

0.003 

1.75 

0.44 

1.367 

6.01E-5 

1.3 

1.510 

0.020 

1.46 

0.008 

1.76 

C . 45 

1.365 

1 .41E-4 

1.536 

1.  510 

0.023 

1.40 

0.008 

1.77 

0.46 

1.359 

2.43E-4 

1.3 

1. 46  0 

0.017 

1 .33 

0.008 

1.79 

0.48 

1.351 

3.11E-4 

2.0 

1.420 

0 .003 

1 .26 

0.003 

1.80 

0.49 

1.  347 

1.C7E-3 

2.25 

1.420 

0.010 

1.22 

0.  009 

1.  81 

0.50 

1 .334 

8 . 5 OE -4 

2.5 

1.  420 

0.  012 

1.13 

0.  009 

1.82 

0.51 

1.309 

2.39E-3 

2.7 

1.400 

0.055 

1.13 

0.013 

1.83 

0 .52 

1 .249 

1.56E-2 

3.0 

1.420 

0.022 

1.16 

0.012 

1.84 

0.54 

1.42  9 

1.97E-1 

3.2 

1.430 

0.008 

1.22 

0.010 

1.86 

0.54 

i . 48i 

6.69E-2 

3.3  923 

1.  430 

0.  007 

1.26 

0.013 

1.87 

0.55 

1.439 

1.51E-2 

3.5 

1.450 

0.005 

1 .2  3 

0.011 

1 .38 

0.56 

1.423 

7.17E-2 

3.75 

1.452 

0.004 

1.27 

0.011 

1.90 

0.57 

1.  398 

2.S0E-3 

4.0 

1.  455 

0.  00  5 

1.26 

0.012 

1. 92 

0.58 

1 .388 

3.69E-3 

4.5 

1 • 46  0 

0.013 

1 .26 

0.014 

1.94 

0.59 

1.377 

9.97E-3 

5.0 

1 .450 

0.012 

1.25 

0.016 

1.97 

0.60 

1.  366 

S.57E-3 

5.5 

i .440 

0.018 

1.22 

0.021 

1.99 

0.61 

1.333 

9.31E-3 

6.0 

1.410 

0.022 

1.15 

0.  037 

2.02 

0.62 

1 .306 

7.96E-2 

6.2 

1.  430 

0.  027 

1.14 

0.039 

2.03 

0.63 

1.421 

6.91E-2 

6.5 

1.460 

0.033 

1.13 

0.042 

2.04 

0.63 

1.374 

2.94E-2 

7.2 

1 .4  00 

0.070 

1 .40 

0.055 

2.06 

0.95 

1.343 

2.49E-2 
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TABLE  10.  REFRACTIVE  INDICES  OF  AEROSOL  PARTICLES  (Continued) 


[Ref. 

12] 

7.9 

1.200 

C.  065 

1.15 

0.  040 

2.12 

0.67 

1.324 

2.79E-2 

8.2 

1 . 01 0 

0.100 

1.13 

0.074 

2.13 

0 .68 

1.324 

3.08E-2 

8.5 

1.300 

0.215 

1.30 

0.090 

2.15 

C.69 

1.33  6 

3.36E-2 

8.7 

2.400 

0.290 

1.40 

0.100 

2.16 

C . 69 

1.33  6 

3.56E-2 

9.0 

2.  560 

C.  370 

1.70 

0.140 

2.17 

0.70 

1.373 

3.65E-2 

9.2 

2.200 

0.420 

1 .72 

0.150 

2.18 

0.70 

1.356 

3.71E-2 

9.5 

1.950 

0.160 

1 .73 

0.162 

2. 19 

0.71 

1.  339 

3.63E-2 

9.8 

1.870 

0.095 

1.74 

0.162 

2.  20 

0.72 

1 .324 

3 . 93E-2 

10.  0 

1.  820 

0.090 

1.75 

0.162 

2.21 

0.72 

1.310 

4.06E-2 

10.591 

1 . 76  0 

0.070 

1 .6  2 

0.120 

2.22 

0.73 

1.271 

5.22E-2 

il.O 

1.720 

0.050 

1.62 

0.105 

2.23 

C . 73 

1.246 

7.31E-2 

11.5 

1.670 

0.047 

1.  59 

0.100 

2.24 

C.74 

1.22  7 

1.05E-1 

12.5 

1.  620 

0.053 

1.51 

0.090 

2.27 

0.75 

1.208 

1.90 E-l 

13.0 

1.620 

0.055 

1 .47 

0.100 

2.23 

0.76 

1.221 

2.23E-1 

14.0 

1 .560 

0.073 

1 .52 

0.085 

2.  31 

0.78 

1.267 

2.71E-1 

14.8 

1.440 

0.100 

1.57 

0.  100 

2.  33 

0.79 

1 .307 

2.92E-1 

15.0 

1.420 

0.200 

1.57 

0.  100 

2.33 

0.79 

1.321 

2.97E-1 

16.4 

1. 750 

0.16  0 

1.60 

0.100 

2.3  6 

0 .31 

1 .407 

3 . 3iE  — 1 

17.2 

2. 08  0 

0.240 

1.63 

0.100 

2.38 

C.32 

1.487 

3.41E-1 

18. 0 

1.980 

0.  180 

1.64 

0.115 

2.4C 

G.83 

1.525 

3 .41 E— 1 

18.  5 

1.  850 

0.  170 

1.64 

0.120 

2.41 

0.33 

1.536 

3.39E-1 

2 0.0 

2.120 

0.220 

1 .63 

0.220 

2.45 

0.35 

1.560 

3.24E-1 

21.3 

2.060 

0.230 

1.77 

0.230 

2.46 

C . 36 

1.  568 

3.18E-1 

22.5 

2.000 

0.240 

1.90 

0.230 

2.  48 

0.87 

1.579 

3.16E-1 

25.  0 

1.  880 

C.  280 

1.97 

0.240 

2.51 

0.39 

1.596 

3.13E-1 

27.9 

1.840 

0.290 

1.89 

0.320 

2.54 

0.91 

1 .612 

3.20E-1 

30  .0 

1.320 

0.300 

1.80 

0 . 420 

2.57 

C.93 

1.614 

3.20E-1 

35.0 

1.920 

0.400 

1.90 

0.500 

2.  63 

0.97 

1.597 

3.16E-1 

40.0 

1.  860 

0.  500 

2.10 

0.600 

2.69 

1.00 

1.  582 

5.61E-1 
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r 


TABLE  11. 

WATER  VAPOR  PRESSURE  AND 
DEW  POINT  TEMPERATURES 

DENSITY  FOR  DIFFERENT 

DEW  POINT 

TEMP. 

VAPOR  PRESSURE 

VAPOR  DENSITY 

(DEG.  KEL 

.) 

( MI  LI  BAR  S) 

(GRAMS/M3) 

233 

0.189 

0.176 

238 

0.314 

0.286 

243 

0.509 

0.453 

248 

0.807 

0.705 

253 

1.254 

1. 074 

258 

1.912 

1.605 

263 

2.863 

2.358 

268 

4.215 

3.407 

273 

6.108 

4.  847 

278 

8.719 

6.797 

23  3 

12.273 

9.399 

288 

17.044 

12.830 

293 

23.3  73 

17.500 

29  8 

31.671 

23.050 

303 

42.43  0 

30.380 

308 

56.236 

39.630 

315 

73.777 

51.190 

TABLE  12.  MONTEREY  BAY  AREA  WEATHER  CONDITIONS  (Refs.  16,  17] 
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TABLE  13.  MOLECULAR  ABSORPTION  AND  SCATTERING  COEFFICIENTS 
FOR  YEARLY  AVERAGES  IN  THE  MONTEREY  BAY  AREA 


FREQUE  NCY 

TRANSITION 

MOL  • ABS. 

MOL.  SCAT. 

TRANS.  ALONG 

ID 

COEF. 

COEF. 

1 KM.  PATH 

( 1 /CM  ) 

( 1/KM) 

(1/KM) 

m 

944.195 

P(20> 

CO  2 

0.1520 

— 

85.  86 

942.3  84 

P ( 22) 

CO  2 

0.1460 

- 

86.44 

2471.340 

P ( 11) 

OF 

0.0230 

- 

97.72 

2496.610 

P(10) 

DF 

0.0184 

- 

98.17 

2527.470 

P ( 12) 

OF 

0.  0172 

- 

98.29 

2553.970 

P C 1 1 ) 

OP 

0.  0199 

- 

99.03 

2580.160 

PCI  0» 

OF 

0.02  93 

- 

97.11 

2605.870 

P ( 9 ) 

DF 

0.0240 

- 

96.65 

2631.090 

P<  8 ) 

DF 

0.0218 

- 

97.34 

2638.396 

PI  11) 

DF 

0.2210 

- 

80.17 

2655. 97C 

P<  7 ) 

DF 

0.0531 

- 

94.83 

2665  .2  OC 

P(10) 

DF 

0.0357 

- 

9 6.49 

2680. 28C 

P (6  ) 

DF 

0. 0547 

- 

94.68 

2703.980 

P ( 5 ) 

DF 

0.02  08 

- 

97.94 

2717.536 

P ( 8 ) 

DF 

0.0930 

- 

91  .12 

2743.028 

P ( 7 ) 

DF 

0.0246 

- 

97.57 

2767.914 

P(  6 ) 

DF 

0.0594 

- 

94.23 

2792.437 

P(  5 ) 

OP 

0. 0439 

- 

95.70 

3373.460 

P( 6 ) 

HF 

0.1240 

- 

88.37 

343  5.1  7C 

P ( 3 ) 

HF 

0.3470 

- 

70.71 

3483. 63  C 

P ( 7 ) 

HF 

2.1800 

- 

11.32 

3531.310 

P ( 6 ) 

HF 

86.5000 

- 

0.00 

3577. 47C 

P ( 5 ) 

HF 

35.0000 

- 

0.00 

3593.800 

P C 8 ) 

HF 

67.2000 

- 

0.00 

3622.710 

P ( 4) 

HF 

35.6000 

- 

0.00 

9433.962 

0.000014 

0.0003283 

99.52 

TABLE  14.  MOLECULAR  ABSORPTION  AND  SCATTERING  COEFFICIENT 
AT  1.06  MICROMETERS  IN  THE  MONTEREY  BAY  AREA 


MOL • ABS.  MOL.  SCAT.  TOTAL  MOL.  TRANS.  ALONG 


COEF. 


COEF. 


EXTINCTION!  1 KM.  PATH 


(1/KM) 

( 1 /KM  ) 

(1/KM) 

(?) 

JANUARY 

9 .2 TOE-6 

9.393E-4 

8.  486E-4 

99.92 

FEBRUARY 

9 .7 00E-6 

8.377E-4 

8.474E-4 

59.92 

MARCH 

9.710E-6 

8.  3 49  E— 4 

8.446E-4 

99.92 

APRIL 

1 0.900E-6 

8. 328E-4 

8.437E-4 

99.92 

MAY 

11 . 800E-6 

8.272E-4 

8.3S0E-4 

59.92 

JUNE 

13.300E-6 

8.  223 E— 4 

8.361E-4 

99.92 

JU  LY 

13.900E-6 

8.232E-4 

8.371E-4 

99.92 

AUGUST 

14.500E-6 

8. 220 E- 4 

8.3  65E-4 

59.92 

SEPTEMBER 

14.100E-6 

8.172E-4 

8.313E-4 

99.92 

OCTOBER 

12.500E-6 

8.236E-4 

8 .3615-4 

99.92 

NOVEMBER 

10.400  E-6 

8.30*E-4 

8.408E-4 

99.92 

DECEMBER 

9 .41 OE-6 

8.344E-4 

8.  428E-4 

59.92 

YEARLY 

11.4005-6 

8.283E-4 

3.397E-4 

99.92 
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TABLE  15.  MOLECULAR  ABSORPTION  COEFFICIENTS  AT  3.8007 
MICROMETERS  IN  THE  MONTEREY  BAY  AREA 


MOL.  ABS. 

TRANS.  ALONG 

ORDER  OF 

COEF  . 

1 KM.  PATH 

BEST  MONTH 

(1/KM) 

(?) 

TO  PROPAGATE 

JANUARY 

0.0194 

98.08 

2 

FEBRUARY 

0.0199 

93.  02 

4 

MARCH 

0.0199 

98.03 

3 

APRIL 

0.021  5 

97.38 

6 

MAY 

0.0223 

97.79 

7 

JUNE 

0.0242 

97.61 

9 

JULY 

0.0249 

97.  54 

11 

AUGUST 

0.0256 

97.48 

12 

SEPTEMBER 

0.0247 

97.  56 

10 

OCTOBER 

0.0227 

97.75 

8 

NOVEMBER 

0.0202 

93.  00 

5 

DECEMBER 

0.  01  93 

98.09 

1 

YEARLY 

0.0218 

98.34 

< 

« 
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TABLE  16.  MOLECULAR  ABSORPTION  COEFFICIENTS  AT  10.591033 
MICROMETERS  IN  THE  MONTEREY  BAY  AREA 


MOL.  ABS. 

TRANS.  ALONG 

ORDER  OF 

COEF. 

1 KM.  PATH 

8EST  MONTH 

(1/KM) 

(?) 

TO  PROPAGATE 

JANUARY 

0.122 

33.51 

z 

FEBRUARY 

0.123 

33.  01 

4 

MARCH 

0.128 

33.02 

3 

APRIL 

0.147 

36.37 

6 

MAY 

0.159 

85.29 

7 

JUNE 

0.135 

33.15 

9 

JULY 

0.195 

32.32 

11 

AUGUST 

0.204 

31.54 

12 

SEPTEMBER 

0.193 

82.43 

10 

OCTOBER 

0.165 

84.78 

8 

NOVEMBER 

0.133 

37.  53 

5 

DECEMBER 

0.121 

98.57 

1 

YEARLY 

0.152 

35.  36 
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TABLE  18.  AEROSOL  EXTINCTION  COEFFICIENTS  FOR  MARITIME 

MODEL  WITH  75%  SEA  SPRAY  AND  25%  RURAL  AEROSOLS 


WAVELENGTH 
(MI  CROMETERS  ) 


AEROSOL  EXTINCTION  AEROSOL  EXTINCTION 

COEFFICIENT  (1/KM)  COEFFICIENT  (1/KM) 

CALCULATED  REPORTED  [Ref.  12] 


2.0 

0.11330 

0.11530 

2.5 

0.09895 

0. 09962 

2.7 

0. 08258 

3.0 

0.10430 

0.10426 

3.2 

0.1080  0 

3.3*523 

0.  10250 

3.5 

0.09945 

0.09899 

3.75 

0 .092S9 

0.09191 

4.0 

0.08802 

0.  08670 

4.5 

0.07875 

5.5 

0 .06007 

0.05928 

6.0 

0. 05555 

0.05485 

6.5 

0.05825 

7.2 

0 .04304 

0.04758 

7.9 

0.  04105 

0.  04063 

8.2 

0.03958 

0.03960 

8.5 

0.03572 

0.  04045 
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FIGURE  4.  RELATIVE  INTENSITIES  OF  HF  LASER  WAVELENGTHS  PRODUCED 


13  13  13 


FIGURE  5.  HF/DF  LASER  SYSTEM  BLOCK  DIAGRAM 


■ 


1 Vacuum  pump 

2 Charcoal  filter 

3 Scrubber 

4 Scrubber-vacuum  valve 

5 By-pass  valve 

6 Cavity  pressure  gauge 

7 Cavity  vacuum  valve 

8 Optical  resonator  mirrors 

9 Brewster  windows 

10  Cavity 

11  Electric  discharge  tube 

12  H2  or  D2  gas  tank 

13  SFg,  02/  He  gas  tanks 

14  Ballast  resistors 

15  Ballast  resistor  cooling  fans 

16  D.C.  power  supply 

17  Water  tank 

18  Cooling  water  pump 

19  Oil  filter 

20  Gas  outlet  to  atmosphere 


FIGURE  5.  HF/DF  LASER  SYSTEM  BLOCK  DIAGRAM  (Continued) 
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X-SCALE-2.Q0E+00  UNITS  INCH. 

T-SCALE=5.00E-02  UNITS  INCH. 

FIGURE  9.  AEROSOL  EXTINCTION  COEFFICIENT  AS  A FUNCTION  OF  WAVELENGTH  FOR 
MARITIME  MODEL  WITH  75%  SEA  SPRAY  AND  25%  RURAL  AEROSOLS 


RELATIVE  HUMID|TV  = 77% 
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X-SCRLE=2.00E+00  UNITS  INCH. 

Y-SCRLE=5.00E-02  UNITS  INCH. 

FIGURE  11.  AEROSOL  EXTINCTION  COEFFICIENT  AS  A FUNCTION  OF  WAVELENGTH 
FOR  DIFFERENT  RELATIVE  HUMIDITIES 


ADD  +2.62E+03  UNITS  TO  ALL  X VALUES. 

X-SCALE=5.00E+00  UNITS  INCH.- 
Y-SCALE=2.00E+01  UNITS  INCH. 

FIGURE  12.  TRANSMITTANCE  ALONG  10  KM  PATH  WITH  HIGH  SPECTRAL  RESOLUTION 
OF  0.01  CM-1  FOR  MONTEREY  YEARLY  AVERAGES  (Near  3.8  Wn) 


ADD  +9.35E+02  UNITS  TO  ALL  X VALUES. 

X-SCALE=5.00E+00  UNITS  INCH. 

Y -SC ALE  = 1 . QOE+O 1 UNITS  INCH. 

FIGURE  14.  TRANSMITTANCE  ALONG  10  KM  PATH  WITH  HIGH  SPECTRAL  RESOLUTION 
OF  0.01  CM"*  FOR  MONTEREY  YEARLY  AVERAGES  (Near  10.6  Mm) 
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=)DD  +2.62E+03  UNITS  TO  ALL  X VALUES. 

X-SCALE-5. OOE+OO  UNITS  INCH. 

T -SC ALE -2 . OOE  +0 1 UNITS  INCH. 

FIGURE  17.  TRANSMITTANCE  ALONG  10  KM  PATH  WITH  DEGRADED  RESOLUTION  FOR 
P = 1018  MB,  T = 290  DEGREES  KELVIN,  PH  ~ = 7.4  8 MB 
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